Draft version March 27, 2013 

Preprint typeset using I^T^^X style emulatcapj v. 12/16/11 



o 



(N 

o 

u 

o 



> 
o 

O 

m 



MODEL-FREE MULTI-PROBE LENSING RECONSTRUCTION OF CLUSTER MASS PROFILES 

Keiichi Umetsu"' 

Draft version March 27, 2013 

ABSTRACT 

Lens magnification by galaxy clusters induces characteristic spatial variations in the number counts of 
background sources, amplifying their observed fluxes and expanding the area of sky, the net effect of 
which, known as magnification bias, depends on the intrinsic faint-end slope of the source luminosity 
function. The bias is strongly negative for red galaxies, dominated by the geometric area distortion, 
whereas it is mildly positive for blue galaxies, enhancing the blue counts toward the cluster center. We 
generalize the Baycsian approach of Umetsu et al. (ApJ, 729, 127, 2011) for reconstructing projected 
cluster mass profiles, by incorporating multiple populations of background sources for magnification 
bias measurements and combining them with complementary lens distortion measurements, effectively 
breaking the mass-sheet degeneracy and improving the statistical precision of cluster mass measure- 
ments. The approach can be further extended to include strong-lensing projected mass estimates, thus 
allowing for non-parametric absolute mass determinations in both the weak and strong regimes. We 
apply this method to our recent CLASH lensing measurements of MACS J1206. 2-0847, and demon- 
strate how combining multi-probe lensing constraints can improve the reconstruction of cluster mass 
profiles. This method will also be useful for a stacked lensing analysis, combining all lensing-related 
effects in the cluster regime, for a definitive determination of the averaged mass profile. 

Keywords: cosmology: observations — dark matter — galaxies: clusters: general — galaxies: clusters: 
MACS J1206. 2-0847 — gravitational lensing: weak — gravitational lensing: strong 



1. INTRODUCTION 

Galaxy clusters, the largest self-gravitating systems 
in the universe, represent a fundamental class of astro- 
physical objects, which contain a wealth of information 
about the initial conditions of primordial density fluctu- 
ations, the emergence and growth of nonlinear structure 
over cosmic time. Clusters can therefore provide funda- 
mental constraints on the nature of dark matter (here- 
after DM) (Clowe et al. 2006), alternative gravity theo- 
ries (Narikawa & Yamamoto 2012), and models of struc- 
ture formation (Allen et al. 2008; Vikhlinin et al. 2009). 

The standard A cold dark matter (ACDM) model 
provides testable predictions for the structure and en- 
vironment of collisionless CDM halos. In this con- 
text, A^-body simulations have established an approx- 
imately self-similar form for the spherically-averaged 
density profile {p{r)) of quasi-equilibrium CDM ha- 
los (Navarro et al. 1997, hereafter Navarro-Frenk- White, 
NEW) over a wide range of halo masses and radii, 
with some intrinsic variance associated with assem- 
bly bias and dynamical structure of individual halos 
(Jing & Suto 2000; Tasitsiomi et al. 2004; Graham et al. 
2006; Fahenbacher et al. 2009; Lapi & Cavaliere 2009; 
Navarro et al. 2010; Gao et al. 2012). 

Massive clusters act as powerful gravitational lenses, 
producing various detectable effects via shifting, mag- 
nifying, and distorting the images of distant back- 
ground sources (Bartelmann & Schneider 2001). Grav- 
itational lensing thus offers a unique opportunity to 
study the underlying matter distribution in and around 
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cluster-sized halos (Umetsu 2010; Kneib & Natarajan 
2011), irrespective of the physical nature, composi- 
tion, and state of lensing matter (Okabe & Umetsu 
2008), providing a direct probe for testing well-defined 
predictions (Oguri & Hamana 2011; Silva et al. 2013). 
Careful lensing work has shown that the total mass 
profiles of clusters exhibit a steepening radial trend 
with a clear profile curvature, in overall agreement 
with the predicted form for the family of CDM- 
dominatcd halos (Gavazzi et al. 2003; Broadhurst et al. 
2005a; Limousin et al. 2007; Umetsu & Broadhurst 2008; 
Lemze et al. 2008; Okabe et al. 2010; Umetsu et al. 
2010; Sereno & Umetsu 2011; Umetsu et al. 2011a,b, 
2012; Coe et al. 2012; Newman et al. 2012b,a). 

The main attraction of cluster gravitational lensing 
is its ability to make a model-independent determina- 
tion of halo mass profiles'^ over a wide range of clus- 
ter radii, when the complementary effects of weak and 
strong lensing are combined together (Broadhurst et al. 
2005a; Umetsu & Broadhurst 2008; Merten et al. 2009; 
Zitrinetal. 2010, 2011b; Umetsu et al. 2010, 2011a,b, 
2012; Oguri et al. 2012). It has been demonstrated 
by our earlier weak-lcnsing work (Broadhurst et al. 
2005a; Umetsu & Broadhurst 2008; Medezinski et al. 
2010) that, without adequate color information, the mea- 
sured distortion signal can be substantially diluted due 
to the contamination by unlensed cluster members, lead- 
ing to biased mass-profile measurements with underes- 
timated profile concentrations, undcrpredicting the ob- 
served Einstein radius from strong lensing. 

^ We remind the reader that model dependence is unavoidable 
to some extent in scientific analysis. In this work we define the 
term "model independent" to refer to those methods without prior 
assumptions about the functional form of the cluster radial profiles 
and mass distributions. 
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Internal consistency of Icnsing measurements can 
be further tested by measuring the complementary 
magnification effects. Gravitational magnification can 
influence the observed surface number density of 
background sources, enhancing their apparent fluxes 
and expanding the area of sky (Broadhurst et al. 
1995; Umetsu & Broadhurst 2008; Hildebrandt et al. 
2009, 2011; Van Waerbeke et al. 2010; Rozo & Schmidt 
2010; Umetsu etal 2011a, 2012; Huff & Graves 2011; 
Ford et al. 2012; Schmidt et al. 2012; Morrison et al. 

2012) . The former effect increases the observable num- 
ber of sources above the limiting flux, whereas the latter 
reduces the effective observing area in the source plane, 
thus decreasing the number of sources per solid angle. 
The net effect is known as magnification bias, and de- 
pends on the intrinsic faint-end slope of the source lumi- 
nosity function. 

In this paper we pursue the utility of all possi- 
ble lensing information available in the cluster regime, 
for making high-accuracy, model-free measurements of 
the cluster mass profile from a joint likelihood analy- 
sis of multi-probe lensing observations. This extends 
the Bayesian approach of Umetsu et al. (2011a) based 
on the unique combination of weak-lensing distortion 
and negative magnification-bias measurements, effec- 
tively breaking degeneracies inherent in a standard weak- 
lensing analysis based on shape information alone (see 
Schneider & Seitz 1995; Schneider et al. 2000). The 
Bayesian method of Umetsu et al. (2011a) has been ex- 
tensively used to reconstruct the projected mass profile in 
a number of clusters: A1689, A1703, A370, C10024+17, 
RXJ1347-11 (Umetsu et al. 2011a), A383 (Zitrin et al. 
2011b), A2261 (Coe et al. 2012), MACS J1206. 2-0847 
(Umetsu et al. 2012), MACS J0416. 1-2403 (Zitrin et al. 

2013) , and MACS J0717.5+3745 (E. Medezinski et al. 
2013, in preparation). In all cases, we find a good agree- 
ment between independent weak- and strong-lensing 
based mass profiles in the region of overlap. 

Here, we shall make a full use of magnification-bias ef- 
fects, by extending source-count measurements into mul- 
tiple populations of background sources, probing a wider 
range of levels of magnification bias. The strong- and 
weak-lensing regimes contribute quite similar logarith- 
mic coverage of the radial profile for a massive cluster 
acting as a super-critical lens (Umetsu et al. 2011a,b, 
2012). In this work, we generalize the Umetsu et al. 
(2011a) method to combine the central strong-lensing 
information with weak-lensing measurements in a joint 
mass-profile analysis. 

The paper is organized as follows. In Section 2 we 
briefiy summarize the basic theory of cluster weak grav- 
itational lensing. In Section 3 we outline our compre- 
hensive methods for obtaining projected cluster mass 
profiles from multi-probe cluster lensing observations. 
In Section 4 we apply our methodology to our re- 
cent CLASH lensing observations of the massive cluster 
MACS J1206.2-0847 (hereafter MACS1206) at 2 0.439 
(Postman et al. 2012; Zitrin et al. 2012; Umetsu et al. 
2012) to derive a joint mass profile solution from the com- 
bined strong-lensing, weak-lensing distortion, positive 
and negative magnification-bias measurements, demon- 
strating how complementary multi-probe cluster lensing 
measurements can improve the reconstruction of clus- 
ter mass profiles. Finally, summary and discussions are 



given in Section 5. 

Throughout this paper, we use the AB magnitude sys- 
tem, and adopt a concordance ACDM cosmology with 
Qm = 0.3, 17a = 0.7, and h = 7?o/(100kms-iMpc-i) 
0.7. Errors represent a confidence level of 68.3% (Icr) 
unless otherwise stated. 

2. CLUSTER WEAK LENSING THEORY 

2.1. Convergence, Shear, and Magnification 

The central quantity in this work is the lensing con- 
vergence K. For cluster lensing, it is expressed as k(6) = 
S](0)/Ecrit, namely the projected mass density 'S{6) in 
units of the critical surface mass density, 

^""=4^^''' (^('')^Dis/Ds, (1) 

where Di, Dg, and Dig are the proper angular diameter 
distances from the observer to the lens (/), the observer 
to the source (s), and the lens to the source, respectively. 
The distance ratio (3{z) represents the geometric strength 
of cluster lensing for a source at redshift z; /3{z) = for 
unlensed objects, z < zi. 

Weak gravitational lensing is responsible for the 
magnification and shearing of the images of back- 
ground sources (Bartelmann et al. 2001; Umetsu 2010; 
Kneib & Natarajan 2011). The image deformation is de- 
scribed by the 2x2 Jacobian matrix A of the lens map- 
ping, 

^ = (1 - k)I - 710-3 - 72Cri, (2) 

with X the identity matrix, a a {a ~ 1,2,3) the Pauli 
matrices, and (71,72) the components of the complex 
gravitational shear 7 71 +172 = |7|e^"'^ with spin-2 ro- 
tational symmetry (e.g., Bartelmann & Schneider 2001). 

The lens magnification is given by the inverse Jacobian 
determinant, ^i{e) = l/\detA{e)\ = 1/|(1 - k)^ - 
In the weak-lensing limit k, \"f \ ^ 1, the shear alone in- 
duces a quadrupole anisotropy of the background images, 
which can be observed from ellipticities of background 
galaxy images (Kaiser et al. 1995). In general, the ob- 
servable quantity for quadrupole weak lensing is not 7 
but the rerfwced gravitational shear g{0) = 'y{0)/[l — K{6)] 
in the subcritical-lensing regime where det^ > (or 1/g* 
in the negative parity region with det^ < and \g\ > 1). 

2.2. Mass-sheet Degeneracy 

The observable distortion field g{6) is invariant under 
the following global linear transformation: 

k(0) ^ Ak(0) + 1 - a, 7(6/) ^ A7(6>) (3) 

with an arbitrary scalar constant A 7^ 0, which 
is known as the mass-sheet or steepness degener- 
acy (see Schneider & Seitz 1995; Saha & Williams 2006; 
Liesenborgs & De Rijcke 2012). This global transforma- 
tion is equivalent to scaling A{9) with A, A{6) XA{9), 
hence leaving the critical curves dctA{9) = in the 
strong-lensing regime invariant. This degeneracy can be 
broken, for example,^ by measuring the magnification 

^ Alternatively, the constant A can be determined such that the 
fc averaged over the outermost cluster region vanishes, if a suffi- 
ciently wide sky coverage is available. Or, one may constrain A 
such that the enclosed mass within a certain aperture is consistent 
with cluster mass estimates from some other observations (e.g., 
Umetsu & Futamase 2000). 
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effects (Section 3.2), because fi transforms as fi{0) 

2.3. Source Populations and Redshift Distributions 

For statistical weak-lensing measurements, we consider 
populations (a = 1,2, ...) of source galaxies with respec- 
tive redshift distribution functions N(^a){z)-^ The mean 
lensing depth for a given population (a) is given as 



d0/3(z)iV(„)(z) 



(4) 



It is useful to introduce the relative lensing strength of 
a source population with respect to a fiducial source at 
infinite redshift as (Bartelmann & Schneider 2001). 

W(a) = /3(q)//3(2 -> oo). (5) 

Then, the convergence and shear for such a fiducial ref- 
erence source are expressed as K(q)(0) — w^ajKoo (6») and 
7(q)(0) = W(a)'-foo{0), respectively. 

3. CLUSTER LENSING METHODOLOGY 

In this section, we generalize the non-parametric 
Bayesian approach of Umetsu ct al. (2011a) for con- 
structing projected mass profiles of individual clusters 
from multi-probe gravitational lensing observations. The 
Bayesian approach enables a full parameter-space extrac- 
tion of model and calibration parameters. It is of par- 
ticular importance to explore the entire parameter space 
and investigate the parameter degeneracies, arising in 
part from the mass-sheet degeneracy. 

3.1. Tangential Distortion 

We construct radial profiles of the tangential distor- 
tion g+ and the 45°-rotatcd component as a function 
of clustocentric radius (see Umetsu & Broadhurst 2008; 
Okabe et al. 2010). In the weak-lensing limit <? « 7, the 
azimuthally-averaged radial distortion profiles satisfy the 
following identity (Kaiser 1995): 



(6) 



with k{9) the azimuthal average of k in the circular 
annulus of radius 9 and k(< 9) the average conver- 
gence within a circular aperture of radius 9, defined as 
K{<9) = {n9')-'J^g,^^^d'9'n{e'). 

With the assumption of quasi-circular symmetry 
in the projected cluster mass distribution k{9) (see 
Umetsu & Broadhurst 2008), the azimuthally-averaged 
tangential distortion is related to k as 



k{< 9) - k{9) Wg[Kooi< 



1 - k{9) 



1 - WgKooiO) 



where Wg is the population-averaged relative lensing 
strength (sec Section 2.3). In the absence of higher or- 
der effects, weak lensing induces only curl-free tangential 

^ In general, we apply different size, magnitude, and color 
cuts in source selection for measuring the lens distortion and 
magnification-bias effects, resulting in different source redshift dis- 
tributions N{z). In contrast to the former effect, the latter does 
not require source galaxies to be spatially resolved, but does require 
a stringent flux limit against incompleteness effects (Umetsu et al. 
2011a). 



distortions (Equation (6). In practice, the presence of x 
modes can thus be used to check for systematic errors. 

3.2. Magnification Bias 

The source number counts for a given magnitude cut- 
off mcut a-i'O modified in the presence of lensing as 
(Broadhurst ct al. 1995) 

N,,i< TOcut) = No{< mcut)M'-''"'(e) = Nobf,{e;s), 

(8) 

where Nq{< mcut) is the unlensed source counts, ap- 
proximated locally as a power-law with slope s = 
dlogj^Q A'o(< mcut)/rf"T,cut > 0. In the weak-lensing 
limit, b^{9;s) « (5s — 2)k{9), so that the signal-to- 

1/2 

noise ratio (S/N) scales as (mcut)|5s(TOcut) ~ 2| 

(Rozo & Schmidt 2010). For a maximally-depleted pop- 
ulation of sources with s — 0, b^{9) = fi^^{6). The 
net magnification effect on the source counts vanishes 
(6^ = 1) when s = 0.4. 

For a mass profile analysis, we calculate the sur- 
face number density 71^(6*) = dN^{9)/dVl of background 
sources as a function of clustocentric radius, by az- 
imuthally averaging N^{6) (for details, see Umetsu et al. 
2011a). The magnification bias is then expressed as 
n^{9) = na^'^-^'^~'^{9) with = dN^{< mcut)/dn and 
jj, the magnification,^ 



(9) 



From Equations (6) and (9), the magnification profile for 
a given source population can be uniquely specified by 
the projected mass density profile Koo(^)- 

3.3. Color- dependent Magnification Bias 

Deep multi-band photometry spanning a wide wave- 
length range allows us to identify distinct populations 
of background galaxies in object color-color space (e.g., 
B — Rc versus i?c — z'), as demonstrated by recent cluster 
weak-lensing work based on Subaru observations (e.g., 
Medezinski et al. 2010, 2011; Umetsu et al. 2010, 2011a, 
2012). Since a given flux limit corresponds to different 
intrinsic luminosities at different source rcdshifts, source- 
count measurements of distinctly-different background 
populations probe different regimes of magnification-bias 
effects. 

The bias is strongly negative for red background galax- 
ies at (z) ~ 1.1 with a fiat faint-end slope s ~ 0.1, re- 
sulting in a net count depletion (Umetsu & Broadhurst 
2008; Umetsu et al. 2011a, 2012), as dominated by the 



geometric area distortion: 6^ « /.t 



-0.75 



For a depleted 



source population with s ^ 0.4, the S/N increases pro- 
gressively as mcut increases. 

On the other hand, the faint blue population of back- 
ground galaxies, lying at (z) ^ 2 (Lilly et al. 2007; 
Medezinski et al. 2010; Umetsu et al. 2010, 2011a, 2012), 

^ Since magnification is nonlinear with respect to K^a , Equation 
(9) is, strictly speaking, only valid for circularly symmetric lenses, 
or applicable in the weak-lensing regime where magnification is 
linearly related to ftoo. In our multi- probe approach, the number 
of data constraints can be sufficiently large to break degeneracies 
in parameter space. In such a case, one can exclude from the 
joint analysis magnification constraints in the nonlinear regime, 
and check for the presence of systematics. 



4 



tends to have a steep intrinsic slope close to the lens- 
ing invariant one (s = 0.4). For such a population, 
source selection can be optimized to maximize the over- 
all S/N (Rozo & Schmidt 2010), by choosing a brighter 
m^ut (i-G., at a lower rtg) corresponding to a steeper count 
slope, say s ~ 0.5, so that the bias is mildly positive and 
a net density enhancement results: 6^ « jjP-^^. 

Hence, combining the distinct blue and red background 
populations probes a wider range of levels of magnifica- 
tion bias (Broadhurst 1995). The relative magnification 
bias, defined as the ratio of the blue to red galaxy counts, 
scales as &^(0;blue)/5^(6>;red) « ^(5))2.5[s(biuc)-s(rcd)]^ 
For the count slopes quoted above, s(blue) 0.5 and 
s(red) = 0.1, so that 6^(blue)/&^(red) « /i. 

In general, combining independent count measure- 
ments from multiple source populations has two major 
advantages to the cluster lensing analysis: First, it im- 
proves the statistical precision of cluster lensing mea- 
surements. Second, it has the further advantage of re- 
duced sensitivity to the intrinsic angular clustering of 
source galaxies, because distinctly different source pop- 
ulations are spatially uncorrelated (Broadhurst 1995). 
In the linear wcak-lensing regime (galaxy-galaxy and 
galaxy- group lensing), the optimal weighting of magni- 
fication signals from different source populations with 
different count slopes has been explored by several au- 
thors (Menard & Bartelmann 2002; Scranton et al. 2005; 
Hildebrandt et al. 2009). 

3.4. Nonlinear Effect on the Source-averaged 
Weak-Lensing Fields 

In general, the weak-lensing effects in the cluster 
regime are nonlinearly related to the underlying lens- 
ing potential. Hence, in general, the averaging operator 
with respect to the source redshift distribution N{z) acts 
nonlinearly on the redshift-dcpendent components in the 
cluster lensing observablcs. Our methodology here can 
be formally generalized to take full account of the non- 
linear effect on the source-averaged lensing profiles, by 
replacing g+{wg) and fep(w^) in Equations (7), (8), and 
(9) as (Seitz & Schneider 1997) 



9+{Wg)- 



dz g+{w[z])Ng{z) 



dz b^{w[z])N^{z) 



dzNg{z) 



dzN^{z) 



(10) 



(11) 



where Ng{z) and N^{z) are the respective redshift 
distribution functions for the shape-distortion and 
magnification-bias measurements. 

In the mildly-nonlinear regime, it is often sufhcient to 
apply a low-order approximation using low-order mo- 
ments of the source-averaged lensing depth, neglecting 
higher-order correction terms. For details, see Appendix 
A. In particular, when the characteristic mean redshift 
{zg) of source galaxies is sufficiently high compared to the 
cluster redshift Zi (see Bartelmann & Schneider 2001), 
we can safely assume that all sources are at the same ef- 
fective redshift, corresponding to the mean lensing depth 
of the source population (e.g., Medezinski et al. 2010, 



2011; Umetsu et al. 2010, 2011a, 2012). 

3.5. Strong Lensing 

The Einstein radius 6^ is a characteristic size-scale 
of strong lensing, tightly related to the cylinder mass 
it encloses, M{< 6'e) (Narayan & Bartelmann 1996; 
Broadhurst & Barkana 2008; Oguri & Blandford 2009; 
Oguri et al. 2012). It describes the area enclosed by the 
tangential critical curve, within which multiple imag- 
ing can occur due to the high surface mass density of 
the lens. For an axisymmetric lens, the average mass 
density within this critical area is equal to Scritj thus 
enabling us to directly estimate the enclosed mass by 
M(< 0e) = 7r(AeE)'Scrit. 

In general, detailed strong-lens modeling with many 
sets of multiple images allows us to determine the criti- 
cal curves with great accuracy, which then provides ro- 
bust estimates of the projected mass enclosed by the 
critical area with an effective Einstein radius 0e (e.g., 
Zitrin et al. 2011b). Accordingly, the enclosed mass pro- 
file 

M(< 9) = Scdt / ^^{e') d^e' = ^(A^^)'S„it,ooKoo(< ( 

J\e'\<9 

(12) 

at the location around 6'e is less sensitive to model- 
ing assumptions and approaches (see Julio et al. 2007; 
Coc ct al. 2010; Umetsu ct al. 2012; Oguri et al. 2012), 
serving as a fundamental observable quantity in the 
strong-lensing regime (Coe et al. 2010, Section 6.1). On 
the other hand, the radial profile slope of S(6') is ill con- 
strained owing to the mass-sheet (steepness) degeneracy 
(see Equation (3)). Since the size of the Einstein ring 
9e{z) grows with the source redshift z, multiple sets of 
multiple images spanning a wide range of source red- 
shifts can be used to construct a reliable inner mass pro- 
file Af (< 9) (e.g., Broadhurst et al. 2005b; Saha & Read 
2009; Zitrin et al. 2009). 

3.6. Bayesian Cluster Mass Profile Reconstruction 
3.6.1. Multi-probe Cluster Lensing Constraints 

We consider all possible lensing information available 
in the cluster regime, namely enclosed mass estimates 
M{< 9) from strong lensing, tangential lens distortion 
g+{d) and magnification bias ?t-^(^) measurements for 
multiple independent populations of background galax- 
ies: 

{"Ma).^}.="i (" = 1. ^coi) (13) 

with A'coi the number of color-selected background sam- 
ples for source-count measurements. We measure the 
lens distortion and count profiles in the subcritical lens- 
ing regime (i.e., outside the critical curves, 9 > 6'e), using 
the same grid of clustocentric annuli {i = 1,2, N^ji). 
Hence, there are a total of iVtot = (1 + A^coOA^wi + -^si 
lensing constraints including A'si central projected mass 
estimates from strong lensing. 

3.6.2. Joint Likelihood Function 

We construct a discretized mass profile from multi- 
probe cluster lensing constraints, extending the Bayesian 
approach by Umetsu et al. (2011a). In the Bayesian 
framework, we sample from the posterior probability 
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density function (PDF) of the underlying signal s given 
the data d, P{s\d). Expectation values of any statis- 
tic of s shall converge to the expectation values of the 
a posteriori marginalized PDF, P{s\d). The covariancc 
matrix C of s is obtained from the resulting posterior 
sample. With the covariance matrix C, we introduce 
an estimator for the signal-to- noise ratio (S/N) for de- 
tection of s, integrated over the radial range considered 
(Umetsu & Broadhurst 2008): 

(S/^f = J2_ s.C-/s,^s'C-'s. (14) 

In our problem, the signal s is a vector containing the 
binned convergence profile {Koo,i}^i with N = iVwi + 
A^si and the average convergence within the innermost 
aperture radius 0min for strong-lensing mass estimates 

'^oOjmin = '*(< ^min): SO that 

^ {/^oo,min; ^oo,z}i— 1 (-^^) 

specified by (N + 1) parameters. 
The Bayes' theorem states that 

P{s\d) cx P{s)P{d\s), (16) 

where C{s) = P{d\s) is the likelihood of the data given 
the model (s), and P{s) is the prior PDF for the model 
parameters. The joint C{s) function for multi-probe clus- 
ter lensing observations is given as a product of their 
separate likelihoods, 

C = CvilCsl = -Cg£^£sl, (17) 

where >Cwi = ^g^f^ and Csi are the likelihood functions 
for weak and strong lensing, respectively, defined as 

ln£^(.)^-lgg["^(")-'-y.-(^'"("))]'(19) 

where (5+, fip(c(), M) are the theoretical predictions for 
the corresponding observations, and all these profiles can 
be uniquely specified by s (see Appendix B). 

The errors a^^i for g^^i due primarily to the variance 
of the intrinsic source ellipticities can be conservatively 
estimated from the data using bootstrap techniques. The 
errors cr^(a),i for '^^(a),* include both contributions from 
Poisson counting errors and contamination due to intrin- 
sic clustering of each source population (Umetsu et al. 
2011a). The strong-lensing mass estimates Mi and errors 
<7M,i can be derived from detailed modeling of multiply- 
lensed images (sec Section 3.5). The inner mass pro- 
file Mi can be measured at several independent aperture 
radii when multiple sets of multiple images are avail- 
able at various source redshifts.® The covariance between 

If there is no strong-lensing constraint available {N^i = 0), 
7c,ui,i represents the average convergence within the inner radial 
boundary of weak-lensing observations. See Umetsu et al. (2011a). 

* Each set of multiply-lcnscd images constrains the mass enclosed 
within their radii, A/(< 9). 



binned k values naturally arises because they are to sat- 
isfy the observed cumulative mass constraints (Coe et al. 
2010). 

3.6.3. Priors 

For each parameter of the model s, we consider a flat 
uninformative prior with a lower bound of s = 0, that is, 
Kcxj.mhi > and Koo,i > 0. Additionally, we account for 
the calibration uncertainty in the observational parame- 
ters, such as the relative lensing strength , the count 
normalization and slope parameters (^0(0)7 S{a)) for each 
color sample (see Umetsu et al. 2011a): 

c= {wg, ri.o(a),S(„)}^^°i, (21) 

giving a set of (3A^coi + 1) calibration parameters to 
marginalize over. 

In practice, the count normalization and slope param- 
eters can be estimated from the counts in cluster out- 
skirts using wide-field imaging data. The mean depth of 
background samples can be either measured from well- 
calibrated photometric redshifts (Umetsu et al. 2012), or 
estimated from deep multi-band photometry, such as the 
30-band COSMOS database (Ilbert et al. 2009). 

3.7. Implementation 

We implement our method using a Markov Chain 
Monte Carlo (MCMC) approach with Metropohs- 
Hastings sampling, by following the prescription outlined 
in Umetsu et al. (2011a). For Bayesian parameter esti- 
mation, we use the location of the marginal maximum a 
posteriori probability (MMAP) for each model parame- 
ter, using the bisection method in conjunction with boot- 
strap techniques (Umetsu et al. 2011a). The covariance 
matrix C for the discrete mass profile s is estimated 
from the MCMC samples. The method has been tested 
against synthetic shear -|-magnification catalogs from sim- 
ulations of analytical NFW lenses performed using the 
public package glafic (Oguri 2010). The results sug- 
gest that both maximum-likelihood (hereafter, ML) and 
Bayesian MMAP solutions produce reliable reconstruc- 
tions with unbiased profile measurements, so that this 
multi-probe lensing method is not sensitive to the choice 
and form of priors, when the shear and magnification are 
combined and hence the mass-sheet degeneracy is fully 
broken. 

4. APPLICATIONS TO CLUSTER LENSING 
OBSERVATIONS; MACS1206 

In this section, we apply our new method to the re- 
cent CLASH observations of MACS1206 presented by 
Umetsu et al. (2012), who performed a joint analysis 
of weak-lensing distortion and negative magnification- 
bias measurements for a reconstruction of the projected 
cluster mass profile, then compared and combined with 
the inner mass profile derived independently from their 
strong-lensing analyses. Here we conduct a joint likeli- 
hood analysis of the full lensing constraints from strong- 
lensing, weak-lensing distortion, positive and negative 
magnification-bias measurements, demonstrating how 
combining multi-probe lensing constraints can improve 
the mass profile reconstruction. 
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4.1. MACS 1206 

MACS 1206 is an X-ray selected CLASH cluster 
(Postman et al. 2012) with a fairly relaxed appearance 
in optical, X-ray, and Sunyaev-Zel'dovich effect (SZE) 
images, as well as in morphology of its brightest clus- 
ter galaxy (BCG) (Ebeling et al. 2009; Gilmour et al. 
2009; Umetsu et al. 2012). No significant offset is ob- 
served between the BCG, X-ray peak, and DM cen- 
ter of mass determined from detailed strong-lens mod- 
eling (Zitrin et al. 2012; Umetsu et al. 2012). For the 
cluster, a good agreement is obtained between the lens- 
ing. X-ray, and SZE mass estimates (Umetsu et al. 2012; 
Rozo et al. 2012), with a virial mass of Mvir = (1.1 ± 
0.2) X IO^^Mq/i^i (Umetsu et al. 2012), indicating the 
hot gas is not far from a state of hydrostatic equilib- 
rium in cluster potential well. On large scales, the clus- 
ter is embedded in elongated rich large-scale structures 
as revealed by the galaxy and weak-lensing mass maps 
(Umetsu et al. 2012). 

The cluster was observed deeply in 16 filters ranging 
from the UV through the optical to the IR on the Hub- 
ble Space Telescope (hereafter, HST) and in five opti- 
cal passbands, BV Rcicz' , with the wide-field Suprime- 
Cam (34' x 27', Miyazaki et al. 2002) on the 8.2 m Sub- 
aru telescope. The majority of Subaru observations 
were taken as part of the Weighing the Giants project 
(von der Linden et al. 2012). We refer the reader to 
Zitrin et al. (2012) and Umetsu et al. (2012) for details 
of the HST (strong lensing) and Subaru (weak lensing) 
observations, respectively. We define the cluster center 
to be the location of the BCG, following Umetsu et al. 
(2012). 

4.2. Weak-lensing Shear and Magnification Constraints 

We derive radial profiles of lens distortion and 
magnification-bias from a reanalysis of Subaru observa- 
tions of Umetsu et al. (2012). We have two color sam- 
ples (iVcoi = 2) of blue and red background populations 
for magnification measurements, and a full background 
sample of blue-|-red galaxies for shape distortion mea- 
surements, as defined by Umetsu et al. (2012, see their 
Figure 3, Tables 3 and 4) using the background-selection 
method of Medezinski et al. (2010). 

For the red counts to measure the negative bias, we 
have 13252 galaxies with a mean depth of /3(red) = 
0.51 ± 0.02, at the magnitude limit of mcut{z') = 
24.6 mag (Umetsu et al. 2012). For this the normaliza- 
tion and slope parameters are estimated as no (red) ~ 
11 A ± 0.3arcmin-2 and s(rcd) = 0.133 ± 0.047 from 
the coverage-corrected source counts in the outer region, 
9 ^ 10'. For the blue counts to measure the positive bias, 
we find that the central density enhancement decreases 
with increasing mcut{z'), vanishing at a deeper cut of 
'mcutiz') 25 mag, where the count slope is close to the 
lensing-invariant one, s(blue) ~ 0.4. Hence, the mag- 
nitude cutoff TOcut(-z') has been chosen to optimize the 
total S/N integrated over all radial bins of the magnifi- 
cation signal (Section 3.2); for this, we find 2740 galaxies 
with /3(blue) = 0.62 ± 0.06, s(blue) = 0.532 ± 0.105, 
and no (blue) = 1.98 ± 0.13arcmin~^ at mcut{z') — 
24.0 mag. Finally, we have 13123 background galaxies 
with /3(back) = 0.54 ± 0.03 for the distortion measure- 
ments (Umetsu et al. 2012). 



We calculate the respective weak-lensing profiles in 
iVwi = 10 discrete radial bins from the cluster center, 
spanning the range 9 = [0.6', 16'] with a constant log- 
arithmic radial spacing Aln6' ~ 0.328. Combining the 
tangential distortion, blue and red count measurements, 
we have a total of 30 constraints from Subaru weak- 
lensing observations. 

In Figure 1 we show the resulting lens-distortion 
(black) , positive (blue) and negative (red) magnification- 
bias measurements of MACS1206 as a function of clus- 
tocentric radius. We find an integrated S/N of 10.2, 
2.9, and 4.7 for the above respective measurements (de- 
fined in analogy to Equation (14)). The total S/N for 
the combined weak-lensing measurements is estimated 
as (S/N)^L = VlO. 22 2.92 + 4.72 ~ 11.6. 

4.3. Strong-lensing Constraints 

The central mass distribution of MACS1206 has been 
tightly constrained by detailed strong-lensing analyses 
based on CLASH _ffS'T imaging and Very Large Telescope 
(VLT) /VIMOS spectroscopic observations (Zitrin et al. 
2012; Umetsu et al. 2012). There are a total of 50 
multiply-lensed images of 13 background sources iden- 
tified for this cluster (Zitrin et al. 2012), spanning a 
wide range of source redshifts, 1 ^ Zg ^ 5.5, spread fairly 
evenly over the central region, 3" ^ 6* ^ 1'. 

For a source at Zg = 2.54, the tangential critical 
curve encloses an area with an effective Einstein ra- 
dius of 6'e = 28" ± 3"; for another lower-redshift sys- 
tem with Zg — 1.03 (Ebeling et al. 2009), the effective 
Einstein radius of the critical area is 6'e = 17" ± 2" 
(see Zitrin et al. 2012; Umetsu et al. 2012),^ correspond- 
ing to model-independent projected mass estimates of 
M{< 17") = 5.80ti;^^ x W^M^h-^ and M{< 28") = 
l-llia25 lO"Af0 /i-^: These are shown to be in broad 
agreement with the aperture mass measurements (Equa- 
tion (12)) obtained from several complementary strong- 
lensing analyses via a variety of modeling methods and 
approaches (for details, see Umetsu et al. 2012, their Fig- 
ure 7). 

In the present analysis, we use the double Einstein- 
radius constraints on the projected total mass M < 9) 
derived at two distinct source redshifts, Zg = 1.03 and 
2.54 {Ns\ = 2). The total S/N for the combined strong- 
lensing constraints is obtained as (S/N)sl — 6.3. 

4.4. Results 

Here we construct discrete mass profiles of MACS1206 
from our multi-probe cluster lensing observations, 
demonstrating how additional strong-lensing and posi- 
tive magnification-bias information can improve the mass 
profile reconstruction. For our full-lensing analysis, we 
have a total of A^tot = (1 + Nco\)N^\ -\- Ngi = 32 con- 
straints from strong-lensing (M), weak-lensing distortion 
{g+), positive and negative magnification-bias measure- 
ments (n^(blue), n^(rcd)), in the range 6'niin = 6'e(zs — 
1.03) = 17" < 9 < 16' = 6i„iax, with A^coi = 2, 
iVwi = 10, A^si = 2, and N ^ N^i + Nsi = 12. The mass- 
profile model s = {Koo,min, Koo^i}fLi then described by 
A^ + 1 = 13 parameters (see Section 3.6.2). Additionally, 

^ Zitrin et al. (2012) used the position and rcdshift of 32 secure 
multiple images of nine systems to constrain their mass model. 
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wc have 7 calibration parameters (c) to marginalize over 
(Section 3.6.3). The projected cumulative mass profile 
M < 6*) is given by integrating the density profile s (see 
also Appendix B) as 

A/(< 9) =7r(L';0,„in)^I]crit,ooKoo,min 

.9 

+2TTDfi:„it,oo / d\n0'9'^K^{d'). (22) 

The results are shown in Figure 2. First, we use only 
the weak-lensing distortion and depletion constraints for 
our joint likelihood analysis (6'min — 0.6', ^max = 16'; 
Ncoi = 1, A^tot = 20, = A^„i = 10), corresponding to 
the weak-lensing analysis of Umetsu et al. (2012). From 
this, we find a consistent mass-profile solution s (green 
circles) with an integrated S/N of 12.1. The ML solution 
has a reduced Xmin 4. 33 for 5 degrees of freedom (dof). 

Next, the mass-profile solution from our fuU-lensing 
constraints is shown in Figure 2 with red squares 
(Xmin/'iof = 13.3/12). It is demonstrated that, 
when the strong-lensing information is included, the 
central weak-lensing bin Koo(< 0.6') is resolved into 
A'si -|- 1 = 3 radial bins, hence improving the de- 
termination of the inner mass profile. The S/N ra- 
tio in the recovered mass profile is 13.6, which is 
consistent within uncertainties with the total S/N 
obtained from the linearly-combined fuU-lensing con- 
straints, V(S/N):^L + (S/N)^L = Vll-62-h6.32 ~ 13.2. 
We find here a minor statistical improvement from 
adding the blue source counts, because the positive 
magnification-bias effect is only marginally detected in 
our analysis (^ 3(t). In Figure 1 we also display the joint 
Bayesian reconstruction of each observed radial profile, 
ensuring consistency and reliability of our lensing analy- 
sis and methods. 

5. SUMMARY AND DISCUSSION 

In this paper we have developed a non-parametric 
Bayesian method for a direct reconstruction of the pro- 
jected cluster mass profile from a multi-probe clus- 
ter lensing analysis (Section 3.6), combining indepen- 
dent strong-lensing, weak-lensing shear and magnifica- 
tion measurements. This multi-probe approach improves 
the accuracy and precision of the cluster lens recon- 
struction, effectively breaking the mass-sheet degener- 
acy (Section 2.2). This work extends our earlier work by 
Umetsu et al. (2011a) to include multiple populations of 
background sources for magnification-bias measurements 
(Sections 3.2 and 3.3) and enclosed total mass estimates 
in the central strong-lensing region (Section 3.5). This 
flexible method applies to both the strong- and weak- 
lensing regimes for full radial coverage (Section 4). 

Magnification bias depends on the intrinsic proper- 
ties of the source luminosity function as well as the 
source distance, in contrast to the purely-geometric 
color- independent distortion effects. The combination 
of count measurements for distinct blue and red back- 
ground populations probes a wider range of levels of mag- 
nification bias (Section 3.3), boosting the significance 
of cluster lensing measurements (Section 4). Further- 
more, combining spatially-uncorrelated distinct popula- 
tions of background galaxies will help reduce the re- 
construction bias due to their intrinsic angular clus- 
tering (Broadhurst et al. 1995; Broadhurst 1995). It is 



also interesting to note that, unlike the shearing effect, 
magnification is in principle sensitive to the sheet-like 
structure (Section 2.2), so that making accurate mag- 
nification measurements is crucial for a robust statisti- 
cal detection of the two-halo term contribution due to 
large-scale structure associated with the central clusters 
(Oguri & Takada 2011; Oguri & Hamana 2011). 

Adding strong-lensing information to weak-lensing is 
needed to provide tighter constraints on the inner den- 
sity profile (Limousin et al. 2007; Umetsu & Broadhurst 
2008; Oguri et al. 2009; Newman et al. 2009; 
Merten et al. 2009; Oguri et al. 2012). The full mass 
profile constraints on individual clusters, derived from 
a joint weak-|-strong lensing analysis, can be further 
stacked together to increase the statistical precision 
of the average mass profile determination. Recently 
Umetsu et al. (2011b) performed a comprehensive 
stacked lensing analysis of four similar mass clusters 
(Umetsu et al. 2011a), by combining weak-lensing 
derived k{9) profiles with their respective inner k(9) 
profiles derived independently from detailed strong-lens 
modeling. In contrast, the new method developed here 
allows us to construct a joint weak-|-strong lensing 
likelihood function of the underlying k signal, by 
explicitly combining weak-lensing measurements with 
strong-lensing enclosed mass estimates M(< 9), a more 
fundamental observable than k{9) in the strong-lensing 
analysis (Coe et al. 2010). Therefore, our new approach 
permits a direct error propagation and thus a more 
accurate derivation of the reconstruction covariance 
matrix C. 

When combining strong and weak lensing constraints 
in different regimes of signal strength and significance, 
it is crucial to account for possible systematic errors in- 
troduced by inherent modeling and prior assumptions, 
especially for strong-lens modeling due to complex, non- 
linear error propagation. Recently, Umetsu et al. (2012) 
introduced a regularization technique to obtain robust, 
conservative error estimates for the binned mass profile 
derived from strong lensing, by calibrating the covari- 
ance matrix to eliminate very small eigenvalues associ- 
ated with large-scale modes where the constraints are 
weak and essentially driven by the prior. It is straight- 
forward to extend our joint likelihood analysis to include 
in Equation (20) the covariance between radial bins. 

An accurate determination of the cluster density 
profile for full radial coverage is crucial for testing 
DM and alternative-gravity paradigms (Newman et al. 
2012b,a; Narikawa & Yamamoto 2012; Geller et al. 2013; 
Silva et al. 2013). Woo & Chiueh (2009) explored in de- 
tail an extremely-light bosonic dark-matter (ELBDM) 
model with a mass of the order 10~^^ cV, as an alterna- 
tive to CDM to account for the perceived lack of small 
galaxies relative to the ACDM model (Klypin et al. 1999; 
Peebles & Nusser 2010). ELBDM with a de-Broglie 
wavelength of astronomical length scales, if it exists, 
may well be in a ground-state Bose-Einstein condensate 
(BEG) and hence described by a coherent wave func- 
tion, behaving effectively as a single scalar field (Hu et al. 
2000). Intriguingly, Woo & Chiueh (2009) showed that, 
ELBDM halos can form steepening density profiles of 
the form similar to the standard CDM, irrespective of 
whether halos form through accretion or merger. How- 
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ever, during a collision between BECs interesting large- 
scale interference occurs which will differ markedly from 
the behavior of cluster galaxies and standard collisionlcss 
CDM. Therefore, it is important to explore this class of 
DM further via more extensive simulations for provid- 
ing testable predictions against detailed cluster lensing 
observations of both relaxed and merging clusters. 

The CLASH survey (Postman et al. 2012) is particu- 
larly designed to generate such useful lensing data, com- 
bining high-resolution 16-band HST imaging with wide- 
field Subaru observations, for a sizable sample of 20 X- 
ray selected {Tx > 5 kcV) relaxed clusters (Zitrin et al. 
2011b; Coe et al. 2012; Umetsu et al. 2012), free of 
lensing-selection bias, and a lensing-selected sample of 
five high-magnification clusters, the majority of which 
are physically-interacting merging systems (Zitrin et al. 
2011a, 2013, E. Medczinski ct al. 2013, in preparation). 
A stacked analysis of the X-ray selected CLASH clusters 
(with a halo bias factor of 6;, ^ 8), combining all lensing- 
related effects in the cluster regime, is highly desirable 
for a definitive determination of the representative mass 
profile out to beyond the virial radius, where the two-halo 
term is expected to be detectable in the averaged k profile 
(Tinker et al. 2010; Oguri & Hamana 2011; Silva et al. 
2013), providing a firm basis for a detailed compari- 
son with the standard ACDM paradigm and a wider 
examination of alternative scenarios (Gao et al. 2012; 
Narikawa & Yamamoto 2012; Silva et al. 2013). 

We thank the anonymous referee for providing valuable 
comments and constructive suggestions. This work was 
made possible in part by the availability of high-quality 
lensing data produced by the CLASH team. I wish to ex- 
press my gratitude to all members of the CLASH team 
who made the data analyzed here possible. I thank 
Elinor Mcdezinski, Mario Nonino, and Alberto Molino, 
for their valuable contributions to the weak-lensing data 
analysis. I am very grateful to Marc Postman for his 
generous support and encouragement. I acknowledge 
Nobuhiro Okabe, Tom Broadhurst, Doron Lemze, Adi 
Zitrin, Dan Coe, Julian Merten, Sherry Suyu, Margaret 
GcUer, and Tzihong Chiuch for fruitful discussions and 
comments. The work is partially supported by the Na- 
tional Science Council of Taiwan under the grant NSC97- 
2112-M-001-020-MY3 and by the Academia Sinica Ca- 
reer Development Award. 

REFERENCES 

Allen, S. W., Rapctti, D. A., Schmidt, R. W., ct al. 2008, 

MNRAS, 383, 879 
Bartclmann, M., King, L. J., & Schneider, P. 2001, A&A, 378, 361 
Bartelmann, M., & Schneider, P. 2001, Phys. Rep., 340, 291 
Broadhurst, T. 1995, ArXiv e-prints, astro-ph/9511150 
Broadhurst, T., Takada, M., Umetsu, K., et al. 2005a, ApJ, 619, 

L143 

Broadhurst, T., Bemtez, N., Coe, D., ct al. 2005b, ApJ, 621, 53 
Broadhurst, T. J., & Barkana, R. 2008, MNRAS, 390, 1647 
Broadhurst, T. J., Taylor, A. N., & Peacock, J. A. 1995, ApJ, 
438, 49 

Clowe, D., Bradac, M., Gonzalez, A. H., et al. 2006, ApJ, 648, 
L109 

Coe, D., Bemtez, N., Broadhurst, T., & Moustakas, L. A. 2010, 

ApJ, 723, 1678 
Coe, D., Umetsu, K., Zitrin, A., et al. 2012, ApJ, 757, 22 
Ebeling, H., Ma, C. J., Kneib, J.-P., et al. 2009, MNRAS, 395, 

1213 



Faltcnbacher, A., Li, C, White, S. D. M., ct al. 2009, Research in 

Astronomy and Astrophysics, 9, 41 
Ford, J., Hildebrandt, H., Van Waerbeke, L., et al. 2012, ApJ, 

754, 143 

Gao, L., Navarro, J. F., Frenk, C. S., ct al. 2012, MNRAS, 425, 
2169 

Gavazzi, R., Fort, B., Mcllier, Y., Pcllo, R., & Dantel-Fort, M. 

2003, A&A, 403, 11 
Gellcr, M. J., Diaferio, A., Rines, K. J., & Serra, A. L. 2013, ApJ, 

764, 58 

Gilmour, R., Best, P., & Almaini, O. 2009, MNRAS, 392, 1509 
Graham, A. W., Merritt, D., Moore, B., Dicmand, J., & Terzic, 

B. 2006, AJ, 132, 2701 
Hildebrandt, H., van Waerbeke, L., & Erben, T. 2009, A&A, 507, 

683 

Hildebrandt, H., Muzzin, A., Erben, T., et al. 2011, ApJ, 733, L30 
Hockstra, H., Franx, M., & Kuijken, K. 2000, ApJ, 532, 88 
Hu, W., Barkana, R., & Gruzinov, A. 2000, Physical Review 

Letters, 85, 1158 
Huff, E. M., & Graves, G. J. 2011, ArXiv c-prints, 1111.1070 
Ilbert, O., Capak, P., Salvato, M., et al. 2009, ApJ, 690, 1236 
Jing, Y. P., & Suto, Y. 2000, ApJ, 529, L69 

Julio, E., Kneib, J.-P., Limousin, M., et al. 2007, New Journal of 

Physics, 9, 447 
Kaiser, N. 1995, ApJ, 439, LI 

Kaiser, N., Squires, G., & Broadhurst, T. 1995, ApJ, 449, 460 
Klypin, A., Kravtsov, A. V., Valenzuela, O., & Prada, F. 1999, 
ApJ, 522, 82 

Kneib, J.-P., & Natarajan, P. 2011, A&A Rev., 19, 47 

Lapi, A., & Cavaliere, A. 2009, ApJ, 692, 174 

Lemze, D., Barkana, R., Broadhurst, T. J., & Rephaeli, Y. 2008, 

MNRAS, 386, 1092 
Liesenborgs, J., & De Rijcke, S. 2012, MNRAS, 425, 1772 
Lilly, S. J., Le Fevre, O., Renzini, A., et al. 2007, ApJS, 172, 70 
Limousin, M., Richard, J., Julio, E., et al. 2007, ApJ, 668, 643 
Mcdezinski, E., Broadhurst, T., Umetsu, K., Bem'tez, N., & 

Taylor, A. 2011, MNRAS, 414, 1840 
Mcdezinski, E., Broadhurst, T., Umetsu, K., ct al. 2010, MNRAS, 

405, 257 

Menard, B., & Bartelmann, M. 2002, A&A, 386, 784 
Merten, J., Cacciato, M., Meneghetti, M., Mignone, C, & 

Bartelmann, M. 2009, A&A, 500, 681 
Miyazaki, S., Komiyama, Y., Sekiguchi, M., et al. 2002, PASJ, 54, 

833 

Morrison, C. B., Scranton, R., Menard, B., ct al. 2012, MNRAS, 
426, 2489 

Narayan, R. , & Bartelmann, M. 1996, ArXiv e-prints, 

astro-ph/9606001 
Narikawa, T., & Yamamoto, K. 2012, JCAP, 5, 16 
Navarro, J. F., Frenk, C. S., & White, S. D. M. 1997, ApJ, 490, 

493 

Navarro, J. F., Ludlow, A., Springel, V., et al. 2010, MNRAS, 
402, 21 

Newman, A. B., Treu, T., Ellis, R. S., & Sand, D. J. 2012a, 

ArXiv e-prints, 1209.1392 
Newman, A. B., Treu, T., Ellis, R. S., et al. 2012b, ArXiv 

e-prints, 1209.1391 
— . 2009, ApJ, 706, 1078 
Oguri, M. 2010, PASJ, 62, 1017 

Oguri, M., Bayliss, M. B., Dahlc, H., et al. 2012, MNRAS, 420, 
3213 

Oguri, M., & Blandford, R. D. 2009, MNRAS, 392, 930 
Oguri, M., & Hamana, T. 2011, MNRAS, 414, 1851 
Oguri, M., & Takada, M. 2011, Phys. Rev. D, 83, 023008 
Oguri, M., Hennawi, J. F., Gladders, M. D., et al. 2009, ApJ, 699, 
1038 

Okabe, N., Takada, M., Umetsu, K., Futamase, T., & Smith, 

G. P. 2010, PASJ, 62, 811 
Okabe, N., & Umetsu, K. 2008, PASJ, 60, 345 
Peebles, P. J. E., & Nusser, A. 2010, Nature, 465, 565 
Postman, M., Coe, D., Bem'tez, N., et al. 2012, ApJS, 199, 25 
Rozo, E., Bartlett, J. G., Evrard, A. E., & Rykoff, E. S. 2012, 

ArXiv e-prints, 1204.6305 
Rozo, E., & Schmidt, F. 2010, arXiv c-prints, 1009.5735 
Saha, P., & Read, J. L 2009, ApJ, 690, 154 
Saha, P., & Williams, L. L. R. 2006, ApJ, 653, 936 



Umetsu 



9 



Schmidt, F., Lcauthaud, A., Massey, R., ct al. 2012, ApJ, 744, 
L22 

Schneider, P., King, L., & Erbcn, T. 2000, A&A, 353, 41 
Schneider, P., & Seitz, C. 1995, A&A, 294, 411 
Scranton, R., Menard, B., Richards, G. T., et al. 2005, ApJ, 633, 
589 

Seitz, C, & Schneider, P. 1997, A&A, 318, 687 
Sereno, M., & Umetsu, K. 2011, MNRAS, 416, 3187 
Silva, L. J. B. e., Lima, M., & Sodrc, L. 2013, ArXiv e-prints, 
1301.1684 

Tasitsiomi, A., Kravtsov, A. V., Gottlober, S., & Klypin, A. A. 

2004, ApJ, 607, 125 
Tinker, J. L., Robertson, B. E., Kravtsov, A. V., et al. 2010, ApJ, 

724, 878 

Umetsu, K. 2010, ArXiv e-prints, 1002.3952 

Umetsu, K., & Broadhurst, T. 2008, ApJ, 684, 177 

Umetsu, K., Broadhurst, T., Zitrin, A., Medczinski, E., & Hsu, L. 

2011a, ApJ, 729, 127 
Umetsu, K., Broadhurst, T., Zitrin, A., et al. 2011b, ApJ, 738, 41 



Umetsu, K., & Futamasc, T. 2000, ApJ, 539, L5 
Umetsu, K., Medezinski, E., Broadhurst, T., et al. 2010, ApJ, 
714, 1470 

Umetsu, K., Medezinski, E., Nonino, M., et al. 2012, ApJ, 755, 56 
Van Waerbeke, L., Hildebrandt, H., Ford, J., & Milkeraitis, M. 

2010, ApJ, 723, L13 
Vikhlinin, A., Kravtsov, A. V., Burenin, R. A., et al. 2009, ApJ, 

692, 1060 

von der Linden, A., Allen, M. T., Applegate, D. E., et al. 2012, 

ArXiv e-prints, 1208.0597 
Woo, T., & Chiueh, T. 2009, ApJ, 697, 850 

Zitrin, A., Broadhurst, T., Barkana, R., Rephaeli, Y., & Bem'tez, 

N. 2011a, MNRAS, 410, 1939 
Zitrin, A., Broadhurst, T., Umetsu, K., et al. 2009, MNRAS, 396, 

1985 

— . 2010, MNRAS, 408, 1916 

Zitrin, A., Broadhurst, T., Coe, D., et al. 2011b, ApJ, 742, 117 
Zitrin, A., Rosati, P., Nonino, M., et al. 2012, ApJ, 749, 97 
Zitrin, A., Meneghetti, M., Umetsu, K., et al. 2013, ApJ, 762, L30 



.0 



-I 



+ 

^ 0.8h 
c 

o 

't^ 0.6h 



0.4 

"D 

to 

c 0.2 
_l 

„ 



I 

c 

"e 

u 
to 

■i-< 

c 

D 
O 

u 



10- 



5- 











R[h \pc] 

500 1000 1500 2000 2500 3000 3500 4000 



T — I — I — r 



T 
« 

# 

# 



T — I — I — r 



T 



T — I — I — r 



T 



T — I — I — r 



T 



T — I — I — r 



T 



T — I — I — r 



T 



1 — I — I — r 



T 



T — I — I — r 



Tangential lens distortion (S/N~10) 
Negative nnagnification bias fronn red counts (S/N-S) 
Positive nnagnification bias from blue counts (S/N~3) j 
Joint solution (shear + magnification + strong lensing) 




I' l l ■ S I 



1 



5 10 
Clustocentric radius, /?[arcmin] 



15 



Figure 1. Cluster weak-lensing radial profiles of MACS1206 {zi = 0.439) obtained from a reanalysis of the Subaru BRcz' data presented 
in Umctsu et al. (2012, see their Figures 4 and 5). The top panel shows the tangential reduced shear profile g+ (squares) based on the 
full baekground sample. The bottom panel shows the coverage-corrected number-count profiles n^, for flux-limited samples of blue and red 
background galaxies (circles). The error bars include contributions from Poisson counting uncertainties and contamination due to intrinsic 
angular clustering of each source population. For the red sample, a strong radial depletion of the source counts is seen toward the cluster 
center due primarily to magnification of the sky area, while a slight enhancement of blue counts is present in the innermost radial bins 
due to the effect of positive magnification bias. Also shown for each observed profile is the joint Bayesian reconstruction from combined 
strong-lensing, weak-lensing tangential distortion, positive and negative magnification-bias measurements (see also Figure 2). 
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Figure 2. Effects of including additional strong-Iensing and positive magnification-bias constraints. The top panel shows the projected 
mass density profile S(f?) of MACS1206 derived from multi-probe cluster lensing measurements. For each case, the innermost central bin 
S(< -Rmin) is marked with a horizontal bar. The green circles show the profile reconstructed from a joint likelihood analysis of 

weak-lensing distortion and negative magnification-bias measurements (Figure 1). The red squares are obtained from a joint analysis of 
the full lensing constraints, including additional strong-lcnsing (boxes, bottom) and positive magnification-bias (Figure 1) measurements. 
The central weak-lensing bin S(< 0.6') is resolved into 3 radial bins when additional double Einstein-radius constraints from strong lensing 
are included. For visual clarity, the reconstructed mass profiles are horizontally shifted by 4% with respect to each other. The red area in 
the bottom panel shows the corresponding enclosed mass profile M < R). The boxes represent the double Einstein-radius constraints on 
the enclosed mass, corresponding to background sources at Zs = 1.03 and 2.54 (Umetsu et al. 2012, see their Figure 7). In each panel, the 
gray area represents the best-fit NFW (-1- mass-sheet) profile derived from a joint analysis of the full- lensing constraints (red squares). 
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APPENDIX 



A. NONLINEAR EFFECT ON THE SOURCE-AVERAGED LENSING FIELDS 

A.L Reduced Shear 

The reduced gravitational shear, g = 7/(1 — k), is nonhncar in k, so that the averaging operator with respect to the 
source redshift acts nonlinearly on k. In general, a wide spread of the redshift distribution of background galaxies, in 
conjunction with the single source-plan approximation, may lead to an overestimation of the gravitational shear in the 
nonlinear regime (Hoekstra et al. 2000). 

Too see this effect, we expand the reduced shear g = g{z) with respect to k = k{z) as 

00 

g = 7/(1 - k) ^ W7oo(l - wKoo)~^ = wyoo (wkoq) . (Al) 

fc=0 

Hence, the reduced shear averaged over the source redshift distribution is expressed as 

00 

{9)^looY.{w'+')nl, (A2) 

fc=0 

where the angular brackets denote an ensemble average over the redshift distribution N{z) of background sources. In 
the weak-lensing limit where /too, |7oo| ^1, (ff) ~ {w)joo = (7)- The next order of approximation is given by 

y,.;^(MUu,-,.l).^-^^^. (A3) 

Scitz & Schneider (1997) found that Equation (A3) yields an excellent approximation in the mildly- nonlinear regime 
of Koo ~ 0.6. Defining /j„ = (w^) / (w)^, a dimensionless quantity of the order of unity, we have the following expression 
for the source-averaged reduced shear valid in the mildly-nonlinear regime: 

with (k) = {w)koo and (7) = {w)'-foD (Seitz & Schneider 1997). For clusters lying at relatively low redshifts, {w^) « (w)^ 
and /„, « 1, leading to the single source-plane approximation: (g) « (7) /(I — ('*))• The level of bias introduced by 
this approximation is A.g/g « {f^ — 1)(k). In typical ground-based deep observations of zi ^ 0.5 clusters (Okabe et al. 
2010; Medezinski et al. 2010; Umetsu et al. 2010, 2011a, 2012), A/^ = /.^ - 1 is found to be of the order of several 
percent, so that the relative error in the g estimate is negligibly small in the mildly- nonlinear regime. 



A. 2. Magnification Bias 

First, let us consider a maximally-depleted sample of background sources with s — dlogiQN{< m)/dm = 0, for 
which the effect of magnification bias is purely geometric, = fJ,~^, and is insensitive to the intrinsic source luminosity 
function. In the nonlinear subcritical regime, the source-averaged magnification bias is expressed as 

{b,) = if^-') = 1 - fl ((7)^ - {n?) - ^.-\{w)) + {fl - 1) {{nf (7)^) , (A5) 

where fi~^{(w)) = (1 — (k))^ — l(7)P- The error associated with the single source-plane approximation is therefore 
^^/i = ifw ~ ^ (7)^) ~ 2A/^((k)^ — (7)^), which is much smaller than unity for background source populations 

of our concern in the mildly-nonlinear subcritical regime ((«) ~ |(7)| ~ O(10~^)). It is therefore reasonable to use the 
single source-plane approximation for calculating the magnification bias of depleted source populations with s ^ 0.4. 

In the positive regime of magnification bias (s > 0.4), on the other hand, interpreting the observed lensing signal 
(i.e., the density enhancement) is more difficult, especially in the nonlinear regime where the fiux amplification factor 
is correspondingly large (say, n ^ 1.5). Therefore, it requires detailed information about the intrinsic source luminosity 
function to apply nonlinear corrections due to the spread of the source redshift distribution. In practice, the distant 
blue population of background galaxies has a distinct, well-defined redshift distribution, which is fairly symmetric 
and peaked at a mean redshift of (zg) ~ 2 (Lilly et al. 2007; Medezinski et al. 2010; Umetsu et al. 2010), so that the 
majority of faint blue galaxies are in the far background of typical cluster lenses, and that the lensing signal has a 
weaker dependence on the source redshift. In such a case, the single source-plane approximation may be justified. 
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B. DISCRETIZED EXPRESSIONS FOR THE CLUSTER LENSING PROFILES 
B.l. Averaged Convergence 

In this appendix, we aim to derive a discrete expression for the mean interior convergence Kqo{< 0) as a func- 
tion of clustocentric radius 9 using the azimuthahy- averaged convergence Koo{9). In the continuous hmit, the mean 
convergence Koo(< 0) interior to radius can be expressed in terms of Hoo{0) as 

T^ooW = ^ ^ rfln0'0'2Koo(0')- (Bl) 

For a given set of {N + 1) annular radii 9i (/ = 1, 2, + 1), defining N radial bands in the range 6'min ^ 9i < 6 < 
0N+1 = ^max, a discretized estimator for Koo(< 0) can be written in the following way: 

9^ 2 

Koo(< ei) = -^Kooi< ^min) + ^ ^ A hi ^.^f^oo (^0 , (B2) 



i=l 



with Aln^i = — 9i)/9i and 9i being the area- weighted center of the ith annulus defined by 9i and 6*^+1; in the 

continuous limit, we have 



d9'9'y{9U, - 9t) = =- ' ' ■ (B3) 
6 tli + t/i+i 

B.2. Lens Distortion and Magnification 

We derive expressions for the binned tangential distortion g^{9;Wg) and magnification ^,{9;w^) in terms of the 
binned convergence Koo, using the following relations: 

g+(6'i;wg) = ^ -=- , (B4) 

H(9^■w^,)^- I — = -, (B5) 

[1 - w^Koo(6'0J 9+{di;Wf,)\ 

where both the quantities depend on the mean convergence Kqo interior to the radius 9i , which is the center of the ith 
radial band of (see Appendix B.l). By assuming a constant density in each radial band and by noting that 

9i is the median radius of the ith radial band, k(< 9i) can be well approximated by (Umetsu et al. 2011a) 

Koo(<^») = ^ Koo(< 6'j) +Koo(< 6l,+ i) , (B6) 

where Koo(< 9i) and Koo(< 0i+i) can be computed using the formulae given in this appendix. 

Accordingly, all relevant cluster Icnsing observables, gj^{9)^nf^{9)^ and M(< 9), can be uniquely specified by the 
discrete convergence profile, s = {Koo.min, too,i}iIi '^ith Koo.min = Koo(^min) and Koo.i = Hoo{9i). 



